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Phase structure of (2 + l)d strongly coupled lattice gauge theories 

Costas G. Strouthos ^ 

^Department of Physics, Box 90305, Duke University, Durham, North Carohna 27708, USA. 

We study the chiral phase transition in (2 + l)d strongly coupled U{N) lattice gauge theories with staggered 
fermions. We show with high precision simulations performed directly in the chiral limit that these models undergo 
a Berezinski-Kosterlitz-Thouless (BKT) transition. We also show that this universality class is unaffected even 
in the large N limit. 
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1. Introduction 

The behavior of symmetries at finite tempera- 
ture is one of the most outstanding problems in 
field theory. In recent years we have witnessed 
revised interest in the chiral phase transition in 
QCD. The problem of symmetry breaking and its 
restoration is intrinsically non-perturbative and 
therefore most of our knowledge about the phe- 
nomenon comes from lattice simulations. How- 
ever, computing quantities in lattice QCD with 
massless quarks is a notoriously difficult problem, 
because most known algorithms break down in 
the chiral limit. In addition, the simulations near 
Tc must be performed on lattices with large spa- 
tial sizes in order to control large finite size effects 
due to the diverging correlation length. 

A useful simplification of QCD occurs in the 
strong coupling limit, which retains much of the 
underlying physics of QCD except for large lat- 
tice artifacts. In this limit chiral symmetry break- 
ing and its restoration at finite temperature have 
been studied using large N and large d expan- 
sions [1]. However, since these approaches are 
based on mean field analysis they cannot help in 
determining the universality of phase transitions. 

Interestingly, lattice QCD with one staggered 
fermion interacting with U{N) gauge fields can 
be mapped into a monomer-dimer system in the 
strong coupling limit [2]. Here, we present nu- 
merical results for the finite temperature criti- 
cal behavior of the (2 -I- l)d model. Our data 
were generated with a recently developed very ef- 
ficient cluster algorithm [3], which allows us to 



perform precision calculations in the chiral limit. 
In agreement with expectations from universality 
and dimensional reduction, we show convincingly 
[4] that the chiral phase trasition belongs to the 
BKT universality class [5]. 

The partition function of the model we study 
here is given by 



ln,b] {ij) 



(1) 



and is discussed in detail in [2,3]. Here = 
0, 1, 2, A'^ refers to the number of monomers on 
the site i, 5^^^ = 0, 1, 2, A^ represents the num- 
ber of dimers on the bond {ij), m is the monomer 
weight, z^jj^ = rjfj are the dimer weights. Note 
that while spatial dimers carry a weight 1, tem- 
poral dimers carry a weight T. The sum is over 
all monomer-dimer configurations [n, h] which are 
constrained such that the sum of the number of 
monomers at each site and the dimers that touch 
the site is always A^ (the number of colors). In 
this work we choose = Ly = L. One can study 
the thermodynamics of the model by working on 
asymmetric lattices with Lt L and allowing T 
to vary continuously. 

2. Results 

In this section we present the results for the 
T 7^ critical behavior of the model. The observ- 
ables used in this work are the chiral susceptibil- 
ity Xc and the winding number susceptibility Xw ■ 
The latter is proportional to the helicity modu- 
lus and describes the response of the system to 
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Table 1 



Fits for Xc and Xw near Tc- 



T 


V 


2r 


Xi/d.o.f 


l/(2c) 


Xi/d.ol 


1.00 


0.222(5) 





0.2 


0.2343(8) 


183.2 


1.02 


0.235(5) 


-0.3(2) 


1.5 


0.2411(5) 


53.4 


1.04 


0.251(5) 


-0.07(2) 


0.4 


0.2483(5) 


2.8 


1.06 


0.249(5) 


-0.03(2) 


0.5 


0.2583(5) 


33.8 


1.10 


0.388(5) 


-0.12(2) 


3.6 


0.2831(5) 


770.0 


1.14 


0.569(6) 


-1.24(3) 


480 







a perturbation that distorts the direction of the 
spontaneous magnetization. The winding num- 
ber susceptibility has been used successfully in to 
demonstrate BKT behavior in other models [6]. 
We fixed = 4 and computed Xc and Xw as a 
function of L. If Tc is the critical temperature, 
then the BKT theory predicts 

Xc oc I L^-'^ [log(L)]0-i25 T = Te (2) 
[ constant T > 

and 

r l/[277(r)] + aiL-"^ T<T^ 
Xw=l [2 + l/log(L/Lo)] T = T, (3) 
[ as cxp (— a4i) T > Tc 

in the large L limit. The critical exponent ??(T) 
is expected to change continuously with T, but 
remains in the range < ri{T) < 0.25. In order 
to confirm these predictions we computed Xc and 
Xw for lattices ranging from L = 32 to L = 750 
and for = 1, ...,32. 

Let us first discuss our results for = 1. In fig- 
ures 1 and 2 we plot Xx; and Xw as functions of L 
for difi^crcnt values of T. Wc find that Xc fits well 
to the form bL^-"^ {\og{L))-^'' when T < 1.06. We 
also find that the logarithmic term is unimpor- 
tant for T < 1.0, whereas at T = 1.04 the value 
2r = 0.7(2) is close to the BKT prediction which 
is 0.125. The values of rj, 2r and the quality of the 
fits Xi/d.o.f arc shown in table 1. We also fit the 
data for Xw to the form {c+ 1 / log{L / Lq)) . The 
values of l/(2c) are also shown in the table. Fi- 
nally, using the fact that x,„ = (2 + 1/ log(L/_Lo)) 
is exactly valid at T(, we fit the data for Xw to 
this form for various values of T. The values of 
xi/d.o.f are shown in the last column of table 
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Figure 1. Plot of Xc vs. L for various T. 



1. Based on where the minimum in xi/d.o.f oc- 
curs we estimate Tc = 1.040(5). Our value of 
77 = 0.250(5) at Tc is in excellent agreement with 
the BKT prediction. 
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Figure 2. Plot of 

Xw vs. L for various T. 
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Figure 3. Plot of 77 as a function of N for 
T = 1.0. The data fit very well to the form rj = 
0.169(6)/Ar + 0.050(9)/Af2 ^i^h xVd.o.f= 1.2. 



Wc also checked that Xc and Xw show similar 
evidence for a BKT transition at larger values of 
N. Using techniques similar to the ones we used 
for A'' = 1 wc computed for various values of 
N. We find that = 0.708(6)iV + 1.40(4) - 
1.07(4)/A'' fits our results very well for all values 
of N with a x^/d.o.f of 1.1. The dependence of 
the coefficients of this polynomial on Lt is still 
under investagation. 
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Figure 4. Plot of as a function of T/Tc for 
various values of N. 



With regards to the N dependence of our re- 
sults wc find two interesting obcrvations. Wittcn 
argued that when the symmetry is U{1) the large 
A'' analysis is still applicable since r] ^ 1/N at 
large A'' [7]. Our results for a fixed T = 1.0, shown 



in figure 3, do agree with his conjecture. Inter- 
estingly, as A'' becomes large and T/T^. is held 
fixed instead of T, wc find that rj ^ Q even in the 
large N limit. Figure 4 shows that r) approaches 
an interesting function of T/Tc as A'' becomes 
large. Extending this observation to QCD, we 
think that the t'Hooft limit (large N with g^N 
held fixed) may be quite similar [8]. 

3. Summary 

We presented high precision results from sim- 
ulations of (2 -I- l)d U {N) strongly coupled lat- 
tice gauge theories at T 7^ and we showed con- 
vincing evidence that the models undergo a BKT 
phase transition. In addition, we showed that this 
universality class is unaffected even in the large 
N limit, implying that the mean field analysis of- 
ten used in this limit breaks down in the critical 
region. 
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